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Abstract 

A wide variety of vacua, and their cosmological realization, may provide an explanation 
for the apparently anthropic choices of some parameters of particle physics and cosmology. 
If the probability on various parameters is weighted by volume, a flat potential for slow- 
roll inflation is also naturally understood, since the flatter the potential the larger the 
volume of the sub-universe. However, such inflationary landscapes have a serious problem, 
predicting an environment that makes it exponentially hard for observers to exist and 
giving an exponentially small probability for a moderate universe like ours. A general 
solution to this problem is proposed, and is illustrated in the context of inflaton decay and 
leptogenesis, leading to an upper bound on the reheating temperature in our sub-universe. 
In a particular scenario of chaotic inflation and non-thermal leptogenesis, predictions 
can be made for the size of CP violating phases, the rate of neutrinoless double beta 
decay and, in the case of theories with gauge-mediated weak scale super symmetry, for the 
fundamental scale of supersymmetry breaking. 



1 Introduction 



Our complicated world is based on hundreds of elements that have different properties. The 
hundreds of elements reflect the existence of hundreds of stable nuclei. The variety of stable 
nuclei requires fine choices of parameters such as c^qcd? «qed ; m u, m d and m e . Even a deviation 
of a few % in q;qcd from the current value could destabilize or stabilize such nuclei as 2 H, 2 He 
and di-neutrons, completely changing the thermonuclear processes, chemical abundances of the 
universe and lifetime of stars 0. Such fine tuning cannot be explained by the naturalness 
principle. Rather, it is easily explained by a combination of the many-universe idea (a cosmo- 
logical diversity in the choice of theories and their parameters) with anthropic selection [21 0] . 
It is only in the sub-universes with fine-tuned parameters that observers made of many species 
of atoms exist. This combination also provides a simple solution to both the why-small and 
why-now problems of the cosmological constant |3J [S] . 

Once we accept that the combination of many-universes with anthropic selection plays a role 
in determining the values of observed parameters, then we need to change our way of thinking. 
The combination not only confirms that the observed values of parameters are consistent with 
the existence of observers, but, because the notion of naturalness is replaced by probability, 
it changes to some extent the questions that should be asked. The fine tuning of parameters 
and initial conditions of inflation is one prominent example: if the probability distribution is 
weighted by the volume of the sub-universes, such a fine tuning is rather probable, as long 
as theory space allows the inflaton mass to scan j^j, or provided appropriate initial conditions 
are realized somewhere in the universe [7j. We call such a space of theories containing lots of 
inflationary regions an inflationary landscape. 

The problem of inflationary landscapes (HI El is that we are no longer able to choose the pa- 
rameters of the inflation model by hand to fit the observed data; rather these parameters (and 
consequently the normalization of the density perturbations) are predicted by a combination 
of statistics from the landscape, cosmological dynamics and anthropic conditions [3]. In par- 
ticular, it is generically true that as the inflaton mass becomes smaller, so the inflated volume 
becomes exponentially larger. This exponential behaviour of the volume-weighted probability 
distribution implies not only that a small inflaton mass is very probable jB], but also that the 
inflaton mass will become as small as possible until it is prevented by some anthropic reason. 
Thus, a naive prediction of inflationary landscapes would be that the inflaton mass is so small 
that the environment is exponentially hard for observers to exist jS] because, for example, the 
density perturbations are too small or too large (this was called "<r problem" in [8]), the baryon 
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asymmetry is too small or there is not enough matter. Whatever the reason, this picture ap- 
parently does not lead to our universe. 1 Since the slow-roll volume factor selects vacua with 
smaller and smaller inflaton mass, we call this the inflation runaway problem. 

The a problem may be solved [H E] if the density perturbations arise from fluctuations of 
light fields other than the inflaton ^TJ Even then, the inflaton mass has to be predicted 
at a moderate value, or otherwise, there could be 2 the same problem as above, for instance, 
in the baryon asymmetry or dark matter abundance. Thus, this is quite a generic problem of 
inflationary landscapes. Note, however, that the probability distribution on observable param- 
eters depends sensitively on the boundary values of unobservable parameters and amplitudes 
being scanned. Thus, the distribution is not exponentially sensitive to observables jT3| in some 
landscapes where the observable-unobservable correlation is sufficiently week jH]. 

In this article, we present an idea of how to stop the runaway behaviour of the inflaton mass, 
along with its concrete implementation. In section l2~Tl we introduce the runaway problem, and in 
section l2~2*l we point out a caveat that allows a solution. It is crucial to distinguish anthropically 
relevant parameters from fundamental parameters, and further 

(a) to have a sharp function, or threshold behaviour, of the former in terms of the latter, and 

(b) to have an anthropic reason that exponentially disfavours some range of relevant param- 
eters. 

The runaway problem can be overcome when both (a) and (b) are satisfied. In section 12.31 we 
see that (a) can be realized as the reheating temperature Tr depends upon the inflaton mass. 
If the inflaton decays dominantly to two heavy particles, Tr drops rapidly as the inflaton mass 
becomes close to twice of the mass of the heavy particle. The inflaton mass is predicted to be 
close to twice the mass of the heavy particle. In this article we assume that only the inflaton 
mass and the cosmological constant are scanned in the landscape, not the mass of the heavy 
particle. 3 Section El explains how Tr can be a relevant parameter. If the baryon asymmetry 
is generated by thermal leptogenesis, a very low reheating temperature leads to a low baryon 
asymmetry and, if the baryon asymmetry is too low, galaxies cannot form before protons decay. 
Thus, too low a Tr is exponentially disfavoured. Therefore, both (a) and (b) can be realized, 
and the runaway problem is solved. In sectionHJwe study the possibility that the heavy particles 

1 After understanding astrophysics and biology better, however, the naive prediction may turn out to be the 
solution to Fermi's paradox "Where is everybody?" 

2 When the inflaton mass does not affect late time cosmology except for the number of e- foldings, as in the 
hybrid inflation model in such a problem does not exist. 

3 We do so as a first attempt to address the runaway problem in the simplest and easiest situation. There 
may be room to lift some of these strong assumptions, but this is beyond the scope of this article. 
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produced by inflaton decay are identified with one of the right-handed neutrinos, hoping for 
constraints and predictions in the neutrino sector. The last section is devoted to conclusions 
and discussion, together with a summary of our predictions. The appendix investigates one of 
the uncertainties in the analysis in sections El and 01 the effect of evaporation of particles from 
proto-galaxies. 

2 The Runaway Problem and How to Solve It 
2.1 The Inflation Runaway Problem 

The probability distribution dV(£)/d£ describes the fraction of observers in the overall universe 
who are in sub-universes where the fundamental parameters, £, are in the range £ to £ + d£ 0. 

dV(0 = dU(ZMOA(a(Z)). (1) 

The initial volume distribution includes the density of vacuum states and the probability 
distribution for the initial amplitudes of the inflaton field. The density of states can be, in 
principle, calculated from a given landscape, or from a top-down theory The volume 

factor V(£) fHl E5] accounts for the expansion of each sub-universe during and after slow-roll 
inflation 

V(0 <xe m °®. (2) 

With these definitions, any volume expansion of the universe prior to the slow-roll stage is 
incorporated in the initial factor /(£). For instance, if a landscape supports eternal inflation, 
its volume expansion is to be incorporated in the initial factor /(£). The anthropic factor A(a) 
depends only on parameters a of the low energy effective theory and of late time cosmology. 
We call these the relevant parameters, some of which are relevant for providing an environment 
suitable for observers. Ultimately they depend on the fundamental parameters, The 
anthropic factor vanishes for many values of the relevant parameters, for example A(m e = 
5 MeV) = 0. However, in this paper the only aspect of the anthropic factor that we consider 
is the formation of galaxies having an appropriate number of baryons, since we pay attention 
only to the scanning of inflation parameters. 

It is rather likely that a landscape supports some eternally inflating regions [T71 E3 E] ; 
our standard-model vacuum with a minute cosmological constant may result from one of them. 
An eternally inflating vacuum can tunnel or jump into regions of the landscape that support 
slow-roll inflation through bubble nucleation or through quantum fluctuation of inflaton fields, 
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eventually ending up in a radiation dominated standard-model universe after reheating (c./. 
|19j). The volumes of eternally inflating vacua keep expanding forever, and hence the volume 
created through nucleation or fluctuation is infinite as well. There is an ongoing effort to 
regularize these infinities so that the probability distribution ((TJ) is well-defined [THl I2D1 EH EE] • 
In this article, we assume that such a regulation is possible. We further assume 4 that, when 
the probability distribution is made well-defined, the classical volume expansion factor (J2J) is 
not exactly cancelled by /(£). Such a cancellation appears unlikely, and one of our motivations 
to consider the scanning of some inflation parameters is to seek an origin for a flat inflaton 
potential in the volume-weighted probability distribution. 

Since the number of e- foldings of slow- roll inflation N e ^ 60 for our sub-universe, V(£) 
ensures that the probability distribution involves a powerful exponential dependence on the 
parameters of the inflation potential After integrating over all fundamental parameters 
and initial amplitudes of various fields that we cannot observe directly, £||, one obtains the 
probability distribution of relevant parameters (fraction of observers who see them) 



Within the anthropic window, where the anthropic factor is changing only mildly, the overall 
probability distribution is governed by the factor in the square bracket. Unless there is an 
accidental cancellation between /(£) and V(£), or unless the region of integration determined 
by a landscape is finely designed, 5 the distribution of observables is generically exponentially 
sensitive to inflation-related quantities such as the density perturbations and the reheating 
temperature [S] . An immediate consequence is that an exponentially small fraction of observers 
in the universe would see values of these parameters in the middle of the anthropic window 
jHJini- Since we do not accept this conclusion that our sub-universe is so extraordinarily special, 
there must be something wrong with these assumptions or arguments. 

In fact, the exponential behaviour from the volume expansion factor V(£) is so steep that the 
first question to be answered is whether the overall probability distribution can be normalized. 
If X(£) and «4.(£) are power-law functions, the probability distribution grows forever and is not 
well-defined. At a minimum, the exponential behaviour must be cut off so that observable 

4 Another assumption that we make in this article is that the volume distribution I(£)V(£)d£ can be treated as 
a continuous distribution over field theory parameters £. This assumption may be translated into a sufficiently 
dense landscape of vacua, or bubble tunneling rates without large hierarchy among them; for more information, 
see 0. 

5 A landscape approximated by an ensemble of chaotic inflation potentials, with a mass-independent cut off 
on the maximum field value, is such a possibility [Sj. Reference |13| presents an explicit realization. 
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parameters have average values that are well-defined. This may be accomplished in three 
different ways. 

If only a limited number of states are available in a landscape after imposing rigid anthropic 
conditions, 6 the probability distribution is always normalizable. On the other hand, the most 
probable vacuum among them does not necessarily satisfy non-rigid anthropic conditions jH] . 

Another possibility is that the initial volume factor provides another exponential distribu- 
tion, so that the combined volume distribution 7(£)V(£)<i£ has two exponentials counteracting 
each other. Various properties of the most probable vacuum will be determined by how the 
two exponential distributions balance. Without studying landscapes of the fundamental theory 
(including eternally inflation regions), it may be hard to understand why the probable vacuum, 
supposed to describe our universe, happens to fall inside the anthropic window. 

The other possibility is that the probability is rendered well-defined by an anthropic factor 
that provides an exponential suppression. Instead of relying on the initial distribution, about 
which we know very little, a special form for the anthropic factor is required, offering the hope 
of observable consequences. We find this case rather attractive and pursue it in this article. 
The problem of this possibility, however, is that it is only outside the anthropic window where 
the anthropic factor begins to decrease and counter the exponential growth in the inflation 
volume. If we assume that both the volume factor J(£)V(£) ~ e 3Ne ^ ± ' ) and the anthropic factor 
A(a(£)) ~ e~ F ( a (™ are governed by some power functions N e (£) and F(a(£)), then the peak 
of the probability distribution lies roughly where 



This implies that the majority of observers in the universe live in sub- universes having an 
environment exponentially hard for observers to exist: A ~ e~ F ( a ^>> ^ e -o(ioo) ^ ^ The pow- 
erful exponential factor V(£) pushes the expectation value of observables outside the anthropic 
window (see Fig. [TJ jSj. 7 We call this the inflation runaway problem. 

6 It should be noted that an anthropic factor from cosmological origins tends not to provide a rigid cut-off 
in the range of parameters [221 IHI - Any conditions that involve density perturbations fall into this category. 
This is because Lagrangian parameters (or sub-universes) determine only the standard deviation a of density 
fluctuations Q, whereas the conditions are imposed on the real fluctuations Q. The anthropic factor A{u) can 
be calculated from A(Q) by assuming a Gaussian distribution: 



Thus, at most the anthropic factor can provide an exponential suppression. 

7 This problem persists even if density perturbations are generated by the curvaton ^T] or modulated reheat- 
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Figure 1: Schematic picture of the inflation runaway problem. From left to right: the anthropic 
factor A, the volume factor V, and the combined probability distribution VA, all on a loga- 
rithmic scale. Normalization of each factor is arbitrary, and hence there is no importance in 
the absolute height in the figure. The region to the right of the vertical axis is the anthropic 
window, and the peak of the probability distribution lies outside the window. 

2.2 An Idea to Solve the Problem 

There is a caveat in the above argument that we exploit in this article. The volume expansion 
factor is a function of the fundamental parameters of inflation models £, whereas the anthropic 
factor depends on the relevant parameters, The caveat is in the map between them: what 

if the map shows a very sharp behaviour? 

Let us assume that a(£) has a region where it is rapidly varying, i.e., \a'\ = \da/d£±\ is 
large, as shown in Fig. [2 Then the factor (da/d£±)~ 1 in (J3]) becomes small as parameters enter 
this region, as illustrated by a green curve in Fig. El and the volume factor V(£) becomes a mild 
function of a (blue curve in Fig. [5J). Although V is exponentially sensitive to £, the dependence 
on a in this region is much weaker. Thus, the factor in the square brackets in (j3J) can form a 
mild peak, as shown in Fig. El If this particular region of a has an overlap with the anthropic 
window, where A(a) depends only mildly on a, the combined probability distribution may have 
a mild peak inside the anthropic window, as shown in the right of Fig. El This explains how 
inflationary landscapes can predict average values of observables inside the anthropic window. 

In this case, the difficulty in (@J) is avoided in the following way. The derivative of the 
distribution function in Q has to be zero somewhere in the anthropic window, forming a peak 

ing |12| mechanism, so that their normalization is predicted to lie within the window, i.e., the a problem is cured 
as in jSJEl- The density perturbations are not the only observables that have to be fixed within moderate values. 
For example, a sub-universe with a reheating temperature sufficiently low to give an exponentially suppressed 
anthropic factor is far from the one we observe. 
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a = T R log 10 (T fl ) 
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f = (m (j> /2mx) log 10 (m^/2m x ) 
Figure 2: A relevant parameter with a sharp dependence on a fundamental parameter 
£. This is interpreted as the dependence of the reheating temperature on the inflaton mass, 
TR^m^) in section EHfl and this figure corresponds to A' ~ 10~ 3 . 



in the distribution. There, 



9 In a 7 d\n£ ± d\nN e - dlnA n 

+ 01 o, + "ti = 0. (5) 



9 In a d In a d In £j_ 9 In a 

where iV e is the number of e-foldings of a sub-universe of maximum probability. If a' is a power 
function of a, the first term is of order unity. If N e and In ^4 ~ — F are simple power functions 
of £j_ and a, respectively, din N e /d\nC,j_ is also of order unity, and the last term is of order 
—F(a). Since we want our own sub-universe to be a highly probable one, N e is approximately 
the number of e-foldings in our sub-universe, which is much larger than one. On the other 
hand, we would like our vacuum to be predicted inside the anthropic window, and hence F(a) 
should be of order unity. Thus, the second term must also be of order unity at the peak of the 
distribution, in spite of the large value of N e . This is possible due to the sharp distribution, as 
long as 

9 In a „- T , nS 
is satisfied. The more the sub- universe inflates, the sharper is the required behaviour of 
2.3 A Mass Threshold in Inflaton Decays 

After inflation the universe must reheat, so that the inflaton must decay. We assume that 
the inflaton mass scans (£ = m^) and that V(m^) (or the combination I(£)V(£)) increases as 
decreases. As the inflaton mass is scanned, mass thresholds in the inflaton decay rate 
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Figure 3: Left: various factors of the probability distribution, namely, A, V and (da/d^)" 1 , 
on a logarithmic scale. The normalization of each factor is arbitrary, and hence there is no 
importance to the absolute height. Right: combined probability distribution dV/da. Note 
that the peak of the distribution is inside the anthropic window. [In this figure, the plateau of 
the volume factor in the threshold region corresponds to N e ~ 10. This is why the peak T^o is 
close to the upper end of the threshold region.] 



automatically provide a sharp behaviour in the reheating temperature Tr\ thus we take 

to be a single relevant parameter, the reheating temperature T^m^). Suppose the inflaton (f) 

decays dominantly to a pair of heavy fermions, <fi — > ipx + ipx, through the operator 

£ = -\(j)i>xipx + h.c, (7) 
and that the universe is reheated by ipx decaying to a fermion ipy and a scalar Y through 

£ = X'^x^yY + h.c. (8) 
For > 2mx, the 2-body decay rate of the inflaton is 

r(^^ x + ^ x )= A 2 ^ ( 9 ) 

107T 

where y//3 describes the mass threshold, and (3 = 1 — {^Imx/m^) 2 decreases sharply as the 
inflaton mass m$ approaches the threshold 2mx from above. The 2-body decay channel closes 
when < 2mj, forcing the inflaton to decay through a 3-body channel 

<j) ip x + ip*x i>x + ipY + Y, (10) 
where ip* x is virtual, with a decay rate 

rw^ + ^ + n^g)^ (ii) 
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TR{m^)\ m<j>KS 2mx ~ \ rp lfll-1/4 



(14) 



for m x < rn^ < 2m X - The decay rate of the inflaton suddenly decreases by roughly (A'/47r) 2 as 
the inflaton crosses the threshold at ~ 2m x . Since the reheating temperature Tr is given 
by ^/rMpi, it drops rapidly 

from T + = min(A, A')^/m Mpi/167r, (12) 

AA' i 

to T_ = — - y/m+Mpi/ 16n, (13) 

as m$ passes through the threshold at 2m x , as shown in Fig. |21 This illustrates the origin of 
the sharp behaviour in a(£): the late-time cosmology parameter a is Tr and the high-energy 
parameter £ is the inflaton mass, m^. Of course the idea in section |2~21 will work only when 
the anthropic factor depends exponentially onn Tr, so that we will be led to consider theories 
where the baryon asymmetry and/or the dark matter relic density depend on Tr. 

Let us examine more explicitly how the averaged value of low-energy observables are deter- 
mined in the example of a(£) = Tj^ra^) above. The threshold behaviour of Tr is approximated 8 
by 

> 2m x (l + A'7(167r)) 
'm+ < 2m x (l - A' 2 /(16tt)) = 

Here, we have assumed that the reheating process above the threshold is governed by the decay 
of the inflaton <p, rather than the decay of ip X - We focus on the part of the threshold with 
> 2m x , and substitute a(£) = T R {m^) = /3 1//4 T + in (0) to give 

T «,o - -|_ (15) 

for a highly probable universe, such as our own. This is an important result. Our mechanism 
requires the reheat temperature of our sub-universe to lie in the upper half 9 of the threshold 
region, so that Tr > y/T + T_ pa a/ max( A, A') / (4tt)T + . However, even for a huge volume factor 

— 1 /2 

V, this requirement is very mild: max(A, A') < 4ii/N e . In the following sections we find that 
a sufficient exponential suppression of the anthropic factor will provide a somewhat stronger 
constraint on max(A, A'). 

8 The approximation of Tr using 3-body decay also exhibits sharp behaviour because the virtual ipx is 
almost on shell as approaches the threshold from below. The definition of (3 is 2[(m^/2rax) — 1] around the 
threshold. 

9 Hcrc, we assume that the volume factor increases as the inflaton mass decreases, as in Fig. |3J There may 
be another peak in the lower half of the threshold region, as in Fig. |3 but it depends on all the power-law 
components of all the factors in the probability distribution whether this peak exists or not, and which peak is 
higher. We just ignore this peak in this article. 
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It is instructive to see in more detail how the probability distribution is maximized in this 
example of a mass threshold in inflaton decays. Assuming there is no significant structure in 
/(m^dm^ around « 2m Xl the initial distribution is approximated by 

Iim^dm^ oc dim^ - 2m x ) oc fpM d\n(T R /T + ) . (16) 

Meanwhile, assuming that the e-fold number N^m^) does not have significant structure around 
~ 2m x , we adopt the Taylor expansion 

N e ( m<f> ) » N e {2m x ) x(l-c ( m0 2 ~^ X )) « #.(1 - c0), (17) 

with an order one coefficient c. Here N e is assumed to increase as m<p decreases (as in the case 
of landscapes of chaotic-inflation regions 0E1IH]), hence the minus sign. The volume expansion 
factor is now approximated by 
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Vh) « e 3iVe e -3ciV e (T«/T + ) 

The overall probability distribution (ignoring the contribution from the anthropic factor in the 
anthropic window) is 

VK)/K)^ » fa.^VblT+f (T^j dln(T R /T + ) (19) 

has a peak around T R /T + ps (l/TVg) 1 / 4 , reproducing (fTK|) . This is how the typical Tr is deter- 
mined from this distribution, and will be the reheating temperature of our sub-universe Tr if 
our vacuum is to be typical. The inflaton mass is predicted in this framework almost at the 
threshold 2mj with a deviation of order {T R ^/T + ) A ~ 1/^e, and the reheating temperature is 
suppressed by (l/iVg) 1 / 4 from the naive expectation T + because the inflaton mass is very close 
to the threshold. 

The above argument shows how the probability distribution can form a peak within the 
anthropic window. We need to further make sure that this is the only peak, or at least, the 
biggest peak in the probability distribution. For <C m x , T^m^) is not a particularly sharp 
function. T R different by a factor of a few would be a result of different by a similar degree, 
which means the number of e-folding different by a similar degree. Increase of N e by a factor of 
2 implies that the volume factor is multiplied by an extra factor of e 3Ne , since e 3NeX2 = e 3Ne e 3Ne . 
Thus, we still need something that counters this exponential growth of probability distribution, 
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and we expect the anthropic factor to do this job; A(a) ~ e~ F ^ with a = Tr now. If both 
N e (£) and F(a(£)) are approximated by simple power functions of £ = m$ below the threshold, 
or Tr < T_, it is sufficient to make sure that 

F(T.) > N e (20) 

and that F{m^) is increasing faster than N e {m^) as decreases, in order to guarantee that 
there is no other peak outside this threshold region. 

In the following sections, we discuss what kind of cosmological scenarios can provide such 
an exponential anthropic factor A(Tr) « e~ F ^ TR \ 

Independently of particular realization of A(Tr), however, one can see the following. It 
follows from the requirement F(T-) <; N e that T_/T + is sufficiently small. This means that 
there is an upper limit on 

T_ max(A, A') 



(21) 



and hence on 

min(A, A') < max(A, A'). (22) 



v/tVIV(167t) 

Thus, the reheating temperature of our universe Tr j0 given by (fH"|) is also bounded from above, 
once N e and are known. Although the upper bound itself depends on such issues as 
particular cosmological scenarios that determine F(Tr), N e , and m^, yet it is generic that 
we have an upper bound on T R for individual cosmological scenario. We see a particularly 
important consequence in the scenario in section [""] 

In this article, we consider the scanning of a single parameter of the inflaton potential that 
causes both and N e to vary. It is beyond the scope of this article to explore the possibility of 
solving the runaway problem for more complicated inflationary landscapes with more than one 
scanning inflationary parameters. We just note here that the scanning of multiple inflationary 
parameters generically leads to exponential behaviour in multiple directions on the space of 
observable parameters, requiring multiple anthropic reasons for exponential suppression in the 
number of observers. 

Our analysis in this article is based on an assumption that only two parameters are scanned 
in a landscape: the cosmological constant and the mass of the inflaton. All other parameters are 
assumed not to be effectively scanned: some parameters may not be scanned in the landscape 
[2*31 124j . or some of the parameters are scanned in the landscape but have already been pinned 
down to the value of our sub- universe via some anthropic argument. The two parameter 
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scanning of (m^, A 4 ) is certainly minimal in considering the runaway problem, and at least the 
first place to start. 



3 An Exponential Anthropic Factor 

The two crucial ingredients for solving the runaway problem are to have a sharp threshold 
behaviour of a relevant parameter as the inflaton mass changes, and to have the anthropic factor 
A decreasing exponentially with this relevant parameter — not just as a power law. Since we 
have seen in the previous section that the reheating temperature exhibits a sharp behaviour as 
the inflaton mass passes a threshold 2mx, it would be sufficient if, in some range of parameters, 
A depends exponentially on Tr. In the following sections, we show how this occurs in some 
cosmological scenarios. 

The basic idea is the following. The baryon symmetry Yr = ub/s decreases in some cosmo- 
logical scenarios if Tr decreases; thermal leptogenesis j2H] is an example — there is insufficient 
production of right-handed neutrinos in the thermal plasma if Tr is less than the mass of the 
lightest one. Although the number of baryons in a given comoving volume determines the max- 
imum number of galaxies available, and hence contributes to the anthropic factor, the baryon 
asymmetry itself is usually only a power-law function of Tr, and is not enough to counter the 
exponential growth of the volume factor. Some other physics consideration is necessary in order 
for Tr to be exponentially relevant to the anthropic factor. 

After density perturbations grow during the era of matter domination, baryons begin to 
self-gravitate to form bound states. In order for these bound states to evolve into such highly 
non-linear structures as galaxies, the systems must release energy. When the number density 
of baryons is too low, the cooling through bremsstrahlung is so inefficient that the time scale 
for the cooling becomes longer than the lifetime of a proton. This observation leads to an 
exponentially suppressed anthropic factor, as shown below in more detail. 

The baryon asymmetry remains constant in thermal leptogenesis as long as Tr > Mi, where 
Mi is the mass of the lightest right-handed neutrino: 



where T Nl / (M^ / M v \) is n Nl /s, assuming out-of-equilibrium decay and ignoring g*, the effective 
statistic degrees of freedom. The CP asymmetry ecp is given by 




(23) 



^cp — 



3 Im^hi^fMi/M,) 



(24) 



lGlT \hi 
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As Tr drops below Mi, however, the baryon asymmetry from thermal leptogenesis is much 
less than (}23|) . because leptogenesis stops before reheating completes. After the end of the 
slow roll era, the inflaton field oscillates and the temperature of the thermal plasma behaves 
as T 4 ~ TrpIH Thermal leptogenesis stops roughly when T ~ Mi, and the subsequent 

entropy production continues until T ~ T R , diluting the lepton asymmetry by a factor of 
(Tr/Mi) 5 Furthermore, the number density of right-handed neutrinos is less than usual 



because 



Thus, the baryon asymmetry decreases as 



Y ^ Y °»(wJ < 26) 

for 100 GeV < Tr < M\. The baryon asymmetry goes down rapidly 10 as Tr decreases, and 
so does the number of Milky- Way type galaxies in a given comoving volume. However, since 
this is only a power law of Tr, this is not enough to produce the desired exponential anthropic 
factor. 

Let us assume that the dark matter density of our universe does not depend on Tr (when 
Tr is in the range, say, above the electroweak scale). Thermal relics of stable particles that 
weigh of order 100 GeV (WIMPs) are a good enough candidate. Under this assumption, the 
epoch of matter-radiation equality remains almost unchanged. On the other hand, a single 
Milky- Way type galaxy has a fixed 11 baryonic mass, of order M ga i ~ (10 11 -10 12 ) M , and 

10 An implicit assumption here is that the branching ratio of the inflaton decay to right-handed neutrinos is 
small enough. In particular, none of the particles ipx, and Y in the model of threshold behaviour in the 
previous section can be identified with the right-handed neutrinos. 

n The total mass of a galaxy is an important parameter of environment for the existence of observers. At least, 
the gravitational potential of a galaxy must be deep enough so that a supernova explosion cannot blow up the 
galaxy. Furthermore, correlations have been found observationally between luminosity (and mass) and metalicity 
of galaxies, and between the metalicity and sizes and orbiting radii of planets |27j . Thus, there will be a band (or 
at least lower bound) in the mass of galaxies suitable for the existence of observers, once all other cosmological 
and particle-physics parameters are fixed. In our universe, it will contain 1O 11 -1O 12 M0, and the lower bound 
of the band will be no less than 10 8 -10 9 M Q — the constraint from the supernova explosion. The band (or the 
lower bound) itself will be different from ours in a sub-universe different from ours. For instance, the supernova 
constraint will be A/ ga i A/g a i,gN. m in.j where Af ga i i sN.min.(Ivir/?'%>) « 10 -3 Mq with the virial temperature at the 
epoch runaway cooling starts. Constraints from the metalicity will require further astrophysical understanding. 
Because we do not know the anthropic band (lower bound) in universes quite different from ours, in this article 
we simply use the estimate M ga i ~ (1O 11 -1O 12 )M0. These enormous astrophysical uncertainties, however, do 
not change one of our main claims that an exponential anthropic factor follows by requiring galaxies to form 
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hence a galaxy has to collect baryons from much wider comoving volume in sub-universes hav- 
ing small baryon asymmetry. The galactic comoving scale enclosing a baryon of order M ga i 
enters the horizon later if the baryon asymmetry is smaller. The horizon entry becomes as 
late as the epoch of matter-radiation equality if Y B were 12 of order ~ (M ga iT 3 ' jra p MV) ~ 
(10~ 15 -10~ 14 )(T eg /0.6 eV) 3 . For Y B less than that, the galactic comoving scale would be pro- 
portional to Y B 1 ^ 3 , the horizon-entry time to Y^ 1 , and the number density of the baryon 
asymmetry [photons] at horizon entry to Y~J [to Y B ]. This number density is diluted by Q 3 as 
the universe continues to expand, until the density perturbation of the galactic comoving scale 
goes non-linear. The density in over-dense regions increases further by of order 100 during 
the virialization process, and the evolution of gravitationally bound systems decouple from the 
expansion of the rest of the universe, 13 giving proto- galaxies with 

p « 100T e 4 g (r B /10- 14 - 5 ) 2 (0.6eV/T eg ) 6 g 3 , (27) 
n, « (T e ,Q) 3 (F B /10- 14 - 5 ) 2 (0.6eV/T e9 ) 6 , (28) 
n B « 1000n 7 y B , (29) 

at the epoch of virialization for Y B < (10~ 15 -10~ 14 )(T e(? /0.6 eV) 3 . Since the galactic comoving 
scale enters the horizon after matter-radiation equality for such low baryon asymmetry, the 
virial temperature of baryons is roughly given by j2H] 

T vir ~ m p Q, (30) 

which is of order 10 keV for Q ~ 10~ 5 . 

For sub-universes with such a low baryon asymmetry, the condition for runaway cooling [2H| 
of the proto-galaxies, 

r c0 oi > (G NP ) 1/2 , (31) 
is not satisfied. Here, r coo i is defined in terms of the energy-loss rate per particle 

dEi oss /dt = —T coo \T v i T . (32) 
before protons decay, providing the anthropic band (lower bound) of galaxy masses does not decrease as fast as 

3/2 

Yg . Our quantitative predictions for the upper bound of the reheating temperature are not greatly affected 
either. Even if there are uncertainties of several orders of magnitude in the anthropic bound on Yb in (|36() or 
(14511 , the limit on Tr is changed very little because the baryon asymmetry depends so sensitively on Tr in (|26|l 
or H44J1 . Sec also section 

12 For n CD Mh 2 ~ 0.113 and Q B h 2 ~ 0, T eq - 0.63 eV. 

13 Thc cosmological constant has to be small enough, so that the universe remains matter dominated until 
this epoch. This only leads to a power-law suppression in the number of states in landscapes, and we ignore 
such a power-law contribution to the probablity distribution. 
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Proto-galaxies can release energy by emitting photons. The contributions to r coo i from the 
two relevant processes, namely bremsstrahlung p + e — > p + e + 7 and Compton scattering 
e~ + 7 — > e~ + 7, are given by [HHI 

a 3 rig 

1-Brems ~ 9 - . =, (33) 



m 2 e ^T vir /m e 
2 7 tt 3 « 2 T 4 



"7 



Comp ~ Q 1 „ o 

810 mf m P 



(34) 



where T 7 ~ is the typical photon energy. To see that (|31|) is not satisfied, first note that 
Compton scattering is more efficient than bremsstrahlung right after virialization; second, the 
rate (J53j) does not satisfy (f3*T|). 

Although [2H] used the condition (J31|) as an anthropic criterion, the cooling condition (J31|) 
not being satisfied only implies 14 that runaway cooling does not happen right after virializa- 
tion. Instead, quasi-static cooling occurs on a time scale of order r~~j, allowing gravitational 
contraction. During quasi-static cooling, the gravitationally bound system maintains virial 
equilibrium, and the Jeans mass remains equal to the total mass of the bound system [25] , 
Hierarchical cooling jSIj does not occur. Thus, sub-structures of galaxies, such as clumps of 
molecular clouds or hydrogen-burning stars, are not formed during quasi-static cooling. As 
the quasi-static cooling proceeds, the system becomes more dense and hotter, 15 and eventually 
the cooling condition (J31|) is satisfied, so that hierarchical cooling starts, and smaller struc- 
tures begin to grow. 16 Thus, we do not find it really safe to claim that the anthropic factor is 
exponentially suppressed whenever the cooling condition (J3T]) is not satisfied. 

The cooling time scale cannot be arbitrary long; 17 galaxies must form before protons decay, 



14 Thc size of the cosmological constant determines the largest structures with order one density perturbations, 
but it is more appropriate to think that the upper bound on the cosmological constant is set by requiring that 
galaxies in the anthropic mass band be formed. Since the landscape picture prefers the largest anthropically 
allowed cosmological constant, it is possible that structures containing baryons of order M ga i (or the lower 
anthropic mass bound for galaxies) are the largest. If so, the destruction of proto-galaxies through mutual 
gravitational interactions before cooling starts is no longer an anthropic problem |29| . See, however, also 
footnotell7land section [5j 

15 At an earlier stage of gravitational contraction, the source of gravitational potential energy is dominated 

1/3 

by dark matter, and thus the temperature of baryons may not increase as in the standard Tg oc p B relation. 
But after a baryonic core with ps ~ Pcbm is formed, as a result of gravitational contraction of baryons, the 
temperature of the baryons (including electrons) begins to increase. 

16 It is the density at this epoch that should be used in evaluating the anthropic conditions in |32| . ensuring 
that the halo not be too dense. See section [3] 

17 Main-sequence stars and planetary systems around them are considered to have formed from cold clouds 
in galactic discs in our universe. Not all the particles are heated up to the virial temperature at the beginning 
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which we assume happens with a lifetime of order r p ~ 10 36 years: 

Coi < r p « 10 36 yrs. (35) 

Since relativistic particles, namely photons, cannot be trapped in a gravitational potential well 
without sufficient scattering interactions, the number density of photons continues to decrease 
because of the expansion of the entire universe. On the other hand, the number density of 
protons remains constant in the virialized systems. Thus, the cooling time scale is set by the 
bremsstrahlung process. 18 Using P%|) - (}3*m) and (}3*3*j) in the anthropic condition (|35[) . we see 
that 

Hr 20 \V2 x nr 5 \ 2 -' , . 

£ 1. (36) 



y b J V Q 

If Yb ^ 10~ 18 , however, the pp pair annihilation in the early universe is incomplete, 19 so 
that p survive with an abundance of order n p /s ~ 10 -18 , which, by charge conservation, is also 
the e + abundance. The bremsstrahlung cooling rate should be evaluated by using n PtP in (JH3J), 
instead of But, since the cross section of e + e~ pair annihilation is larger than that of (e,p) 
bremsstrahlung, pair annihilation starts before cooling. As the number densities of electrons 
and positrons decrease, the (e,p) bremsstrahlung cooling rate also decreases; the total energy 
carried by electrons (and positrons) is only n e ^/n p , p . In the end, only electrons are left with the 
abundance n e = n p — n p , and Tei-cms reduces to the one with % again. Anti-protons annihilate 
with protons before proto-galaxies begin to cool. Thus, even for Yb ^ 10~ 18 , the anthropic 
bound (jUD is still valid. 20 

The anthropic bound (j^BI has to be satisfied by the actual density perturbation Q. Assuming 
that Q follows a Gaussian distribution characterized by its standard deviation a, we obtain the 
anthropic factor 

A(* t Y B )K [ dQe -(QM^ e - { io--VYB)^(2xio^n\ (37) 

i(Q/2xlO- 5 )>(10-20/y fl )6/5 



of virialization, as argued in jSHJ- But, if it takes too long time for proto-galaxies to cool, would-be discs may 
be heated and evaporate before runaway cooling starts. This anthropic condition is not taken into acount in 
section and 0] See discussion in section 

18 Hcre, we assume that all the particles in proto-galaxies remain bounded gravitationally forever. However, 
this is not really correct; some particles evaporate from proto-galaxies before runaway cooling starts. See 
discussion in the appendix for the effects of evaporation. 

19 TW thanks G. Perez for bringing this issue to our attention. 

20 We assume that e + -e _ and p-p pair annihilation emit photons and relativistic pions that freely escape from 
the gravitational potential of proto-galaxies. 
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Since Yb oc for Tr < M\ in the thermal leptogenesis scenario, this anthropic factor indeed 
depends exponentially on Tjj(m^). 

In order to see whether this exponentially small anthropic factor is small enough to contain 
the exponentially large volume factor, we need to be a little more specific about the inflationary 
landscape. We approximate the inflationary landscape by an ensemble of chaotic inflation 
potentials, with only one parameter to be scanned. The volume factor is given by V ~ e 3JVe 
with N e (m^) ~ Mpi/rritf,, because the largest value of the inflaton field for classical slow-roll 

3 /2 1 /2 

inflation is <p ~ M pl jm^ . The amount of inflation increases as decreases, and thus the 
volume factor can, in principle, be countered by the threshold behaviour and the anthropic 
factor we have discussed. The asymptotic behaviour of the anthropic factor A ~ e~ F ^ TR(jri4 '^ is 
given by 

F{T R {m*)) oc V B 12/ V 2 oc T R ^m^ oc m" 27 ' 2 , (38) 

where a ~ m^/M p i for chaotic inflation and Tr oc a/T(0) oc m 3 J 2 for m x ^ rn^ <C 2m x [see 
(JUJ)]. Thus, as rriff, decreases, F(Tji(m^)) increases much faster than A^(m^) and we conclude 
that the exponential anthropic factor from thermal leptogenesis is able to tame the volume 
factor from chaotic inflation. 

The anthropic factor (|37|) should begin to decrease within the range of XL <Tr < T + covered 
by the threshold behaviour, so that the peak of the probability distribution is determined by 
the argument in the previous section. The requirement ([20)1 on T_ now reads 

Yb (T- ) < 1(T 20 x (m^/M pl ) 5 / 12 « 1(T 22 , (39) 

where the COBE normalization was used to determine the inflaton mass ~ 10 13 GeV 
in chaotic inflation. This is translated into T_/Mi ^ 10 -1 ' 6 , or equivalently, (T_/T + ) ^ 
10- 1 - 6 7V e " 1/4 (AX 1 /T jR>0 ) £ 10^ 3 (M 1 /T /?>0 ), with N e ~ [2.4 x 10 18 GeV/10 13 GeV]. Thus, the 
coupling constants for the reheating have to satisfy only a mild condition: max(A, A') ^ 
10~ 2 (M 1 /T R ). Therefore, we obtain the upper bound on the reheating temperature as outlined 
in the previous section: 

T Rfi < VlO 12 GeVMi, (40) 

and the lightest right-handed neutrino has to be somewhat lighter than 10 12 GeV. 

We have seen in this section that thermal leptogenesis can solve the runaway problem of 
chaotic inflation, with only a mild constraint on the coupling constants involved; A, A' ^ 10~ 2 . 
We have also obtained some predictions; the upper bound on the lightest right-handed neutrino 
mass Mi < 10 12 GeV, and the upper bound on the reheating temperature of our universe (|40p. 
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Although these predictions are not directly testable in low-energy experiments, they provide 
extra constraints on the parameter space of thermal leptogenesis. 

4 Non-thermal Leptogenesis 

The key ingredient in section 2 is the existence of a heavy particle ipx, with mass mx ~ 10 13 
GeV that is not scanned in a landscape. This hypothetical particle ipx realized the sharp 
threshold behaviour in the reheating temperature, and stopped the inflaton mass runaway at 
2mx- Can ipx be related to other observable physics? 

The tiny masses of left-handed neutrinos combined with the see-saw mechanism suggest 
that right-handed neutrinos exist, with masses of order 10 15 GeV for order 1 Dirac Yukawa 
couplings and 10 13 GeV for order 10" 1 Dirac Yukawa couplings. Thus, a natural question 
arises: can a right-handed neutrino be identified with the hypothetical particle ipx*? If this idea 
works, one may be able to make further arguments in the future why mx is not scanned in 
landscapes by looking at how the masses of right-handed neutrinos arise. In this section, we 
show that one of the right-handed neutrinos can indeed be identified with the threshold particle. 
This means that the inflaton mass, and consequently the density perturbation, is determined 
by the physics of right-handed neutrinos. 

Let us assume that the inflaton decays to a pair of right-handed neutrinos iVi with a branch- 
ing ratio very close to 1 via the interaction 

C = -A0AyVi + h.c (41) 

The right-handed neutrinos decay to a lepton / and a Higgs scalar h u through 

£ = Aj a JV i Z a /i tt + h.c, (42) 

These two operators are the same as (J7J) and (jSJ), allowing exactly the same behaviour as in 
section 01 The inflaton mass is predicted to be ~ 2m^ 1) and the reheating temperature of our 
sub- universe is suppressed by roughly relative to the reheating temperature naturally 

expected when the inflaton mass is safely above the threshold 2m^ x . 

An important difference from the scenario in the previous section, however, is that the 
baryon asymmetry does not decrease as fast as in (|26|). The lepton asymmetry from the 
thermal leptogenesis certainly goes down as in (|2(j[). but there is an extra contribution from 
the the leptogenesis in the decay of right-handed neutrinos produced by inflaton decay (22] • Its 
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lepton asymmetry is given by 

Y B ~ 0.5^e CP = Y Bfi ( J*-") , (43) 
^0 \Tr,oJ 

and decreases very slowly as Tr decreases. Thus, the lepton (or baryon) asymmetry is domi- 
nated by this contribution for low reheating temperatures. Since the baryon asymmetry still 
depends on Tr, the sharp behaviour of T^m^) can be used to solve the runaway problem. 
The lepton asymmetry, however, depends so weakly on Tr, as opposed to (|26|). that a baryon 
asymmetry as low as Yb ~ 10~ 22 is not achieved for Tr larger than the electroweak scale. 
For reheating temperatures lower than the electroweak scale, the baryon asymmetry decreases 
much faster, i.e., 

~ AOOGeVW T R \ 5 

because the sphaleron process stops after the temperature of the thermal plasma drops below 
100 GeV, and the baryon asymmetry generated at high temperature is diluted by radiation 
produced from inflaton decay. Thus, the reheating temperature below the threshold, T_, is 
likely to be lower than the electroweak scale, but not much less. 

The anthropic bound (|36|) on Yr (and on the density perturbation Q) was based on an 
assumption that the dark matter energy density does not change, but, with a reheating tem- 
perature lower than the electroweak scale, this assumption may not hold. If cold dark matter is 
a thermal relic of a stable particle with mass of order 100 GeV (WIMP dark matter), production 
after reheating may be highly suppressed. The relic number density is frozen when the plasma 
temperature is around 10 GeV, and is diluted by a factor of order (Tr/10 GeV) 5 due to the 
following entropy production from inflaton decay. Gravitinos produced from thermal scattering 
are also a candidate for dark matter in gauge-mediated supersymmetry breaking scenarios, and 
their number density also decreases as the reheating temperature decreases — not exactly in the 
same manner for Tr above the electroweak scale, though. 

For Tr sufficiently less than 10 GeV for WIMP dark matter, the dominant component of 
dark matter is left-handed neutrinos. The effect of a lower dark matter density is incorporated 
by taking T eq to be m„ ps 0.03 eV in (J2*7)l - (j2*9*|l . Thus, the anthropic bound (|3l)j) is replaced by 



10" 21 \ 3 /2 x 10~ 5X 25 



Yb J V Q 



< 1. (45) 



Note that the galactic comoving scale enclosing a baryonic mass of M ga \ is longer than the 
free-streaming length of neutrinos for Yb as low as 10 -21 , and the neutrinos can be treated as 
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cold dark matter in the earlier stage of galaxy formation. The same argument as in the previous 
section leads to an exponential cut off for low reheating temperature and small inflaton mass, 
and it is sufficient to take T_ so that 

F B (T_) < 1(T 21 x (m^/M pl f 12 « 10" 24 . (46) 

Thus, T_ ~ 1 GeV is sufficient for T Rfi ~ 10 5 5 GeV, and T_ ~ 10 GeV for T Rfi ~ 10 10 ' 5 GeV. 
Since Y R depends on T R quite sensitively in ()44|) for low T R , the exponential anthropic factor in 
f!37|) is powerful enough to tame the volume factor of chaotic inflation. At the same time, one 
will notice that various uncertainties already mentioned in footnotes in the previous section do 
not lead to a large uncertainty in the constraint on T_. 

One can also see that the bound ()45|) (and (|46|)) is valid for gravitino dark matter as well, 
essentially because left-handed neutrinos become the dominant cold dark matter and the other 
components are small and no longer relevant. 

Now we can derive the upper bound on the reheating temperature of our universe, using the 
constraint (|46[). Higher reheating temperature T R0 (and higher T + ) requires larger coupling 
constant, because T + ps min(A, \' i=l a )xl0 15 GeV, where ~ 10 13 GeV is used for chaotic 
inflation. This means larger T_/T + ~ max(A, \')/(4ir). Since the constraint (|4l)Jl sets an 
upper limit on T_, T Rj0 cannot be arbitrary high. Let us assume iV e 1/4 ~ 10. For T R . o ~ 10 
GeV, T + ~ 10 8 GeV, which means that min(A, X' i=1 a ) ~ 10~ 7 . The threshold behaviour 
of T^m^) can cover the range from T + to T_ ~ (max(A, A')/47r)T + <: 1 GeV. Thus, for 
A ~ \' i=1 a ~ 10~ 7 , the anthropic factor is already sufficiently exponentially small at the lower 
end of the threshold. For T R q ~ 10 8 GeV, T_ cannot be less than 100 GeV, and the anthropic 
factor is not exponentially small. The upper bound on the reheating temperature of our sub- 
universe, Tr,o, is only slightly above 10 7 GeV. 

At the same time we obtain the upper bound on the Dirac Yukawa coupling constants of 
the lightest right-handed neutrinos: 21 |A^ =la | ^ 10" 7 . 

Leptogenesis through inflaton decay accounts for the baryon asymmetry of our universe 
Y B ,o ~ (0.4-0.9) x 10~ 10 only when T Rfi > 10 6 GeV, since the CP asymmetry in (J2U) cannot be 
more than order one. Thus, we see that the range of reheating temperatures of our universe is 
narrowed down to one order of magnitude, 10 6 GeV ~ 10 7 GeV. This means that the CP phase 
in leptogenesis is no less than of order 0.1. If there is no significant accidental cancellation with 
other combinations of CP phases in the lepton sector, there is a chance that CP violation will 

21 Dirac Yukawa coupling constants of other right-handed neutrinos are not constrained. Thus, the CP asym- 
metry can be large enough to account for the lepton/baryon asymmetry of order 1CP 10 . 
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be seen in the next-generation of neutrino oscillation experiments. 

The other prediction, | A- =1 Q | ^ 10 , implies that the mass matrix of left-handed neutrinos 
is virtually of rank 2, if there are 3 right-handed neutrinos. Some predictions on the rate of 
neutrinoless double beta decay are available in the rank 2 case; see [33] for more information. 

These predictions are based on chaotic inflation, where the inflaton mass is determined by 
the COBE normalization m$ ~ 10 13 GeV. We immediately learn that one of the right-handed 
neutrinos has a mass 

With weak-scale supersymmetry breaking 22 arising from gauge mediation, our prediction 
for the range of reheating temperatures Tr has a significant further implication. In these 
theories, the gravitino is the best dark matter candidate for a wide range of gravitino masses, 
100 keV-100 GeV. The relic density of gravitinos produced from the thermal plasma is given 

by mm 

n j2 f T R \ AO MeV\ /W100 GeV)\ 2 

and depends on two totally unknown parameters m 3 / 2 and Tr. Only the ratio (TR/m 3 / 2 ) can be 
determined from the observed value of Qh 2 ~ 0.113. But, now that the reheating temperature 
is essentially predicted the gravitino mass is determined, 7723/2 — 20 MeV-200 MeV. This in 
turn determines the fundamental scale of supersymmetry breaking, \fF « 10 8 -10 9 GeV. 



5 Conclusions and Discussion 

If the fundamental theory allows more than one vacuum, and if many vacua are realized cos- 
mologically, then cosmology may play an important role in choosing the low energy theory 
and its parameters. In inflationary landscapes, where models and parameters of inflation are 
scanned cosmologically, sub-universes with a small inflaton mass dominate the volume of the 
universe, making slow-roll inflation a probable consequence. The problem of inflationary land- 
scapes (El EI] , however, is that the probability distribution tends to be exponentially sensitive 
to some late-time cosmological parameters, and hence does not correctly predict their observed 
values. We call this the inflation runaway problem. 

Landscapes and many-universes are about the world outside the horizon of our universe, 
and we usually think that no sign of them can be seen. One also usually considers that the an- 

22 In supcrsymmctric theories, the Affleck-Dine mechanism is another way to generate the baryon asymmetry. 
However, a negative mass squared is generated only when the Kahlcr potential satisfies certain conditions. If a 
landscape supports a region where the conditions are satisfied, an extra consideration is necessary. 
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thropic conditions depend on only a limited number of relevant parameters, most of which have 
already been measured, so that there are no predictions. Without a well-understood top-down 
tool to figure out the statistic distribution of vacua on a landscape, it seems almost impossible 
to determine theories or parameters that are yet to be tested in experiments. It is important 
to note, however, that inflationary landscapes have a generic prediction: the probability dis- 
tribution is exponentially sensitive to some parameters. This is where the runaway problem 
comes from. The problem is so severe that there may be very few possible solutions, from which 
indications may be obtained for particle physics to be explored in the future. 

Some ideas have already been proposed to evade this problem. For instance, eternal infla- 
tion of false-vacuum type or chaotic-type may strongly select only one inflationary region in a 
landscape [HI E] . The observational consequences for the curvature of the universe and for lower 
multipoles of CMB spectrum are discussed in ^J] and references therein. Maybe the funda- 
mental theory of phsycis has an inflationary landscape carefully designed to avoid the volume 
factor exponentially sensitive to some parameters [HI EI] • If such landscapes involve chaotic 
type inflation, we may expect tensor perturbations. A fundamental cut off scale parametri- 
cally larger than the observed Planck scale may also have something to do with the apparent 
fine tuning problems of certain low-energy supersymmetric theories [HH]- An alternative is to 
have density perturbations generated by a curvaton or by modulated reheating [HI E] , whose 
observational consequences are discussed in the literature [TT) IT3]. 

In this article, we have proposed a new solution to the runaway problem, which is connected 
to the particle physics of inflaton decay and baryogenesis. In spite of the exponential distribution 
from the inflationary landscapes, observables are predicted successfully in the middle of the 
allowed anthropic windows. The key ingredient is to distinguish between parameters describing 
inflation and those that are directly observed. The latter are not necessarily given by simple 
power-law functions of the former. For instance, the reheating temperature drops sharply as 
the inflaton mass approaches the threshold for two body decay. In this case, the inflationary 
landscape predicts the inflaton mass to be close to twice the mass of the threshold particle. 

A sufficiently sharp threshold behaviour for Tr requires small couplings in operators relevant 
to reheating, and hence a low Tr in our universe that depends on the cosmological scenario: 
Tr ^ 10 12 GeV in the thermal leptogenesis scenario of section [21 and 10 6 GeV ^ T R ^ 
10 7 GeV in the inflaton-decay leptogenesis scenario of section Chaotic inflation is assumed. 
In both sections, the anthropic limit was set by requiring that bremsstrahlung and runaway 
hierarchical cooling of proto-galaxies start before protons decay. Although we have ignored 
various uncertainties, the upper bound obtained on Tr illustrates a completely new way to 
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obtain information on particle cosmology from anthropic reasoning. 

In section 0] the particle produced by inflaton decay was identified with a right-handed 
neutrino, allowing the inflaton mass, and some other parameters including the normalization of 
density perturbations, to be set by physics of the neutrino sector. This provides a mechanism 
for understanding the values of inflation parameters in terms of observable particle physics. 
In this particular scenario, there are several predictions that again demonstrate how this line 
of anthropic argument can lead to observable tests. One CP violating phase of the neutrino 
sector must be of order 0.1 or larger. The low energy neutrino mass matrix is forced to be 
essentially rank 2, so that the rate for neutrinoless double beta decay will be constrained by a 
particular relationship involving the neutrino mixing angles and the observed neutrino masses. 
Furthermore, in supersymmetric theories with gauge mediation, dark matter is predicted to 
be gravitinos with mass in the range 20-200 MeV. This corresponds to a fundamental scale of 
supersymmetry breaking of 10 8 -10 9 GeV, which can be tested by precise measurements of the 
superpartner spectrum and by the nature of the decay of the next-to-lightest superpartner [37] . 

It should be noted that the above predictions are based on an assumption that the anthropic 
limit is set by the comparison between the rate of bremsstrahlung cooling and proton decay. 
When various other anthropic conditions are imposed, some of which are listed in footnotes in 
sectional the anthropic limit in the Yb-Q plane may be stronger, and the upper bounds on Tr 
may become weaker. Including other uncertainties in galaxy formation and in the anthropic 
conditions, the above predictions may be modified. There is no doubt that further studies in 
astrophysics will be important in making the predictions more precise and reliable. As a step 
in this direction, we discuss possible effects of particle evaporation from proto-galaxies in the 
appendix, concluding that the upper bound on T R will not be greatly changed, partly because, 
in scenarios in section |3] and EJ Yg is highly sensitive to Tr. 

As the astrophysics and anthropic conditions become better understood, the upper bounds 
of the reheating temperature will change numerically. But, the idea for stopping the runaway 
by a combination of threshold effects and exponential anthropic factors will survive; extra 
anthropic conditions would only strengthen our conclusion that the exponential distribution 
of inflationary landscapes can be tamed. No matter how sophisticated and complicated the 
astrophysical analysis becomes, the upper bound on the reheating temeperature persists as 
long as the threshold behaviour of the inflaton decay plays a crucial role. 
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Appendix: Evaporation 

In proto-galaxies that cannot cool, some fraction of the baryons in the high-energy tail of the 
Boltzmann distribution have sufficient speed to escape, as noted in j^UJ. We ignored this effect 
in the main text, and devote this appendix to investigating the extent to which the conclusions 
could be changed by evaporation. 

There are two important points in discussing evaporation from proto-galaxies. One is that 
the evaporation of particles from isolated proto-galaxies results in their gravitational contraction 
|41j . The other is that the proto-galaxies are comprized of multiple components: dark matter 
particles and baryons (including electrons and possibly their anti-particles). We begin for 
simplicity with the evolution of proto-galaxies made purely of cold dark matter particles, and 
later extend the discussion to take account of baryons. 

Nearly 1% of particles in the Maxwell-Boltzmann distribution have velocities larger than 
the escape velocity of proto-galaxies and hence evaporate. The time scale of evaporation is 
set by how quickly the lost population of particles in the high-energy tail of the Boltzmann 
distribution is restored by particle-particle interactions. Two-body interactions via gravity 
always exist, no matter what the nature of the dark matter particle, with an efficiency given 

r ~ < G ^)' / ' 2 (48) 

^ DM 

where N^m is the typical number of dark matter particles in one proto-galaxy Note that 
the number of baryons is very large, of order Nb ~ 10 68 for M ga i ~ 10 n M Q , and Ndm is 
extremely large if dark matter is composed of elementary particles with a moderate mass, such 
as 100 GeV, or even lighter. Thus, the gravitational encounter of two dark-matter particles is 
not effective within the time scale of bremsstrahlung cooling. 
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When dark-matter particles have two-body non-gravitational scattering with a cross section 
(Tod, the relaxation due to this scattering occurs at the rate ubm^doVvw, where udm is the 
number density of dark-matter particles. The evaporation is roughly 1% of the relaxation rate, 
so that evaporation of dark-matter particle does not occur before bremsstrahlung cooling, if 

lO" 2 n D MO-DDV viT ^ r BrC ms- (49) 
Using (|27)l and (jHHJ), this holds true as long as 

for dark matter with Pdm/ s — 0.8 eV, and 

for neutrino dark matter with n v /s ~ 10 _L5 . It is clear that cross sections for axion-axion 
scattering and neutrino-neutrino scattering satisfy (J50|) and ([51)1. respectively, for the values of 
Yb given in the right-hand sides. Thus, the anthropic bounds (|36|) and (|45|) in section |3] and 0] 
are unaffected by the evaporation of dark matter in such cases. 

While the inequality (|5U|) is satisfied by most dark matter candidates, for example WIMPs, 
it could be violated when the interactions between dark matter particles are significant. In 
that case, the anthropic bound (|36|) in section |3] will be affected; dark-matter evaporation 
leads to the gravitational contraction of proto-galaxies, and then the baryon number density is 
also enhanced and bremsstrahlung cooling begins to work. At some point in the gravitational 
contraction due to evaporation and quasi-static bremsstrahlung cooling, runaway cooling starts. 
Thus, evaporation changes the estimate of the time scale before the hierarchical cooling starts, 
which is to be compared with the proton lifetime in (J35|) . However, even for such dark matter 
candidates, the estimated anthropic bound in (J36j) is not greatly affected by the evaporation of 
dark matter particles, since unitarity implies that([50|) is not violated by a large amount — not 
by of order 10 10 . When the anthropic bounds on Yb is translated to an upper bound on the 
reheating temperature Tr, the effect will be quite minor since Yb is very sensitive to Tr in (|26j). 

Let us now bring baryons into the discussion, ignoring the evaporation of dark matter par- 
ticles. Baryons-baryon scattering has a large cross section and the Boltzmann distribution is 
repopulated much more quickly. We do not go into a detailed estimate of the evaporation rate 
of baryons, but an important point is that this evaporation is accompanied by a gravitational 
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contraction of the baryons. As long as there is no efficient energy transfer between dark-matter 
particles and baryons, baryons fall into the potential well of proto-galaxies as the evaporation 
proceeds, so that the evaporation soon stops. Thus, the discussion in the main text is unaf- 
fected by the evaporation of baryons (and electrons) if the energy transfer between dark matter 
particles and baryons does not occur effeciently. 

Consider the case that the dark-matter particle is lighter than the proton, for example 
the case of axions or neutrinos. Since the axion-nucleon and neutrino-nucleon scattering cross 
sections are so small, an inequality similar to (|51|) is safely satisfied. Thus, the energy transfer 
does not take place between axion or neutrino dark matter with baryons before bremsstrahlung 
cooling starts, and the conclusion of the previous paragraph holds. 

When dark matter particle has a mass of order 100 GeV-1 TeV, the kinetic energy of a 
proton typically increases by of order unity in a single collision with a dark-matter particle. 
Once a proton is kicked by that amount, it has a good chance to escape the gravitational 
potential. The typical rate for this to happen is ubm^dnVviv, which is to be compared with 
TBrcms- Energy transfer from dark matter to baryons does not occur before bremsstrahlung 
cooling, if the dark-matter nucleon scattering cross section a^N satisfies 

Most candidates for dark matter have a scattering cross section with protons that is much 
smaller than 0.1 pb, and hence the discussion in sectionElis not affected by baryon evaporation 
from proto-galaxies. 

In the rare cases that energy is transferred from dark matter to baryons, the baryons may 
evaporate efficiently from proto-galaxies since dark matter particles keep supplying energy 
to the baryons allowing them to escape. Since protons and electrons are light and minor 
components in proto-galaxies, the energy supply continues until proto-galaxies loose virtually 
all their protons and electrons. Thus, bremsstrahlung cooling should work in a time scale 
shorter than that of the energy transfer to baryons and evaporation of baryons. This anthropic 
condition sets a constraint in the Y^-Q plane in addition to the proton lifetime constraint (|35|). 
An argument can be constructed essentially in the same way as in the main text, and in the 
end, the main qualitative result is maintained even with evaporation: namely, an exponential 
anthropic factor results from the failure of structure formation, and an upper bound on the 
reheating temperature is required to make Yb(T_) small enough. Since ()52|) is not a stringent 
constraint, and since Yb is very sensitive to Tr, the upper bound on Tr is not affected very 
much, even quantitatively. 
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It should be noted, however, that the discussion so far, both in the appendix and hence in 
the main text, assumes that proto-galaxies are homogeneous, allowing a crude estimate of the 
impact of evaporation. The limit on the cross section, such as 10 pb, may change a little when 
one takes into account the concentration of particles at the centre of the proto-galaxies and 
the profiles of the matter distribution and gravitational potential. The limit may also change 
when additional anthropic conditions, such as those in footnotes in section |31 are taken into 
account, or when M gal ~ 10 11 M Q is replaced by an anthropic lower bound on M gal , yet to be 
determined precisely. Further study of such issues may change the argument in this appendix 
quantitatively, but the way of thinking will remain valid. 

References 

[1] e.g., J. Barrow and F. Tipler, "The Anthropic Cosmological Principle" Oxford University 
Press, 1986. 

[2] E. Baum, "Zero Cosmological Constant From Minimum Action," Phys. Lett. B 133 (1983) 
185; 

S. W. Hawking, "The Cosmological Constant Is Probably Zero," Phys. Lett. B 134, 403 
(1984); 

S. R. Coleman, "Black Holes As Red Herrings: Topological Fluctuations And The Loss Of 
Quantum Coherence," Nucl. Phys. B 307, 867 (1988); 

S. B. Giddings and A. Strominger, "Loss Of Incoherence And Determination Of Coupling 
Constants In Quantum Gravity," Nucl. Phys. B 307, 854 (1988). 

[3] J. Garcia-Bellido, A. D. Linde and D. A. Linde, "Fluctuations of the gravitational 
constant in the inflationary Brans-Dicke cosmology," Phys. Rev. D 50, 730 (1994) 
|arXiv:astro-ph/9312039| ; 

A. Vilenkin, "Predictions from quantum cosmology," Phys. Rev. Lett. 74, 846 (1995) 
|gr-qc/9406010| ; 

J. Garcia-Bellido and A. D. Linde, "Stationarity of inflation and predictions of quantum 
cosmology," Phys. Rev. D 51, 429 (1995) |arXiv:hep-th/9408023| . 

[4] S. Weinberg, "Anthropic Bound On The Cosmological Constant," Phys. Rev. Lett. 59, 
2607 (1987); 

G. Efstathiou, Mon. Not. R. astron. Soc. 274, L73 (1995); 



27 



H. Martel, P. R. Shapiro and S. Weinberg, "Likely Values of the Cosmological Constant," 
Astrophys. J. 492, 29 (1998) |arXiv:astro-ph/9701099| . 

[5] J. Garriga, M. Livio and A. Vilenkin, "The cosmological constant and the time of its 
dominance," Phys. Rev. D 61, 023503 (2000) |arXiv:astro-ph/9906210| ; 
S. A. Bludman, "Vacuum energy: If not now, then when?," Nucl. Phys. A 663, 865 (2000) 
|arXiv: astro-ph / 99071681 . 

[6] A. Vilenkin, in 0. 

[7] A. D. Linde, "Chaotic Inflation," Phys. Lett. B 129, 177 (1983). 

[8] B. Feldstein, L. J. Hall and T. Watari, "Density perturbations and the cosmo- 
logical constant from inflationary landscapes," Phys. Rev. D 72, 123506 (2005) 
|arXiv:hep-th/050623"5] . 

[9] J. Garriga and A. Vilenkin, "Anthropic prediction for Lambda and the Q catastrophe," 
|arXiv:hep-th/0508005l 

[10] S.Webb, u Ifthe Universe is Teeming with Aliens... Where Is Everybody ?' Springer (2002). 

[11] A. D. Linde and V. Mukhanov, "Nongaussian isocurvature perturbations from inflation," 
Phys. Rev. D 56, 535 (1997) |arXiv:astro-ph/9610219| ; 

K. Enqvist and M. S. Sloth, "Adiabatic CMB perturbations in pre big bang string cosmol- 
ogy," Nucl. Phys. B 626, 395 (2002) |arXiv:hep-ph/0109214| ; 

D. H. Lyth and D. Wands, "Generating the curvature perturbation without an inflaton," 
Phys. Lett. B 524, 5 (2002) |arXiv:hep-ph/0110002| ; 

T. Moroi and T. Takahashi, "Effects of cosmological moduli fields on cosmic mi- 
crowave background," Phys. Lett. B 522, 215 (2001) [Erratum-ibid. B 539, 303 (2002)] 
|arXiv:hep-ph/0110096 . 



[12] G. Dvali, A. Gruzinov and M. Zaldarriaga, "A new mechanism for generating density 
perturbations from inflation," Phys. Rev. D 69, 023505 (2004) |arXiv: astro-ph / 030359 1 



L. Kofman, "Probing string theory with modulated cosmological fluctuations," arXiv: 
|astro-ph/03036T4| 



[13] A. Linde and V. Mukhanov, "The curvaton web," arXiv:astro-ph/0511736 



[14] R. Bousso and J. Polchinski, "Quantization of four-form fluxes and dynamical neutraliza- 
tion of the cosmological constant," JHEP 0006, 006 (2000) |arXiv:hep-th/0004134| ; 
M. R. Douglas, "The statistics of string / M theory vacua," JHEP 0305, 046 (2003) 
arXiv:hep-th/0303194| ; 



28 



S. Ashok and M. R. Douglas, "Counting flux vacua," JHEP 0401, 060 (2004) 
|arXiv:hep-th/0307049| ; 



T. Banks, M. Dine and E. Gorbatov, "Is there a string theory landscape?," JHEP 0408, 
058 (2004) |arXiv:hep-th/0309170| ; 

F. Denef and M. R. Douglas, "Distributions of flux vacua," JHEP 04 05 , 072 (2004) 



|arXiv:hep-th/0404116 



A. Giryavets, S. Kachru and P. K. Tripathy, "On the taxonomy of flux vacua," JHEP 
0408, 002 (2004) |arXiv:hep-th/0404243|; 



O. DeWolfe, A. Giryavets, S. Kachru and W. Taylor, "Enumerating flux vacua with en- 
hanced symmetries," JHEP 0502, 037 (2005) ]arXiv:hep-th/0411061| . 

[15] A. D. Linde, "Eternal Chaotic Inflation," Mod. Phys. Lett. A 1, 81 (1986), "Eternally 
Existing Selfreproducing Chaotic Inflationary Universe," Phys. Lett. B 175, 395 (1986). 

[16] A. D. Linde and A. Mezhlumian, "Stationary universe," Phys. Lett. B 307, 25 (1993) 
|arXiv:gr-qc/9304015] ; 

A. D. Linde, D. A. Linde and A. Mezhlumian, "From the Big Bang theory to the theory 
of a stationary universe," Phys. Rev. D 49, 1783 (1994) |arXiv:gr^q c / 9306035 1 . 

[17] P. J. Steinhardt, in The Very Early Universe, Proceedings of the Nuffield Workshop, 
Cambridge, 1982, eds: G. W. Gibbons, S. W. Hawking, and S. T. C. Siklos (Cambridge 
University Press, 1983), pp.251; 

A. Vilenkin, "The Birth Of Inflationary Universes," Phys. Rev. D 27, 2848 (1983). 
[18] L. Susskind, "The anthropic landscape of string theory," arXiv:hep-th/0302219. 
[19] B. Freivogel, M. Kleban, M. R. Martinez and L. Susskind, "Observational consequences of 



a landscape," arXiv:hep-th/0505232 



[20] A. Vilenkin, "Unambiguous probabilities in an eternally inflating universe," Phys. Rev. 
Lett. 81, 5501 (1998) |arXiv:hep-th/9806185| ; 

V. Vanchurin, A. Vilenkin and S. Winitzki, "Predictability crisis in inflationary cosmology 
and its resolution," Phys. Rev. D 61, 083507 (2000) |arXiv:gr-qc/9905097| ; 
J. Garriga and A. Vilenkin, "A prescription for probabilities in eternal inflation," Phys. 
Rev. D 64, 023507 (2001) arXiv:gr-qc/0102090| ; 



M. Tegmark, "What does inflation really predict?," arXiv:astro-ph/0410281 

A. H. Guth, "Inflation and eternal inflation," Phys. Rept. 333, 555 (2000) 

|arXiv: astro-ph/0002 156] . 



29 



[21] J. Garriga, D. Schwartz-Perlov, A. Vilenkin and S. Winitzki, "Probabilities in the infla- 



tionary multiverse," arXiv :hep-th/0509184 



R. Easther, E. A. Lim and M. R. Martin, "Counting pockets with world lines in eternal 
inflation," arXiv: astro-ph /05 1 1 233 

[22] H. Martel, P. R. Shapiro and S. Weinberg in g]. 

[23] N. Arkani-Hamed, S. Dimopoulos and S. Kachru, "Predictive landscapes and new physics 
at a TeV," |arXiv:hep-th/0501082| 

[24] B. S. Acharya, F. Denef and R. Valandro, "Statistics of M theory vacua," 
|arXiv:hep-th/0502060| 

[25] for a review on leptogenesis, see W. Buchmuller, R. D. Peccei and T. Yanagida, "Lepto- 
genesis as the origin of matter," |arXiv:hep-ph"70 502169 , 

[26] E. Kolb, M. Turner, "The Early Universe" Addison-Wesley, 1990. 

[27] N. C. Santos, G. Israelian, M. Mayor, R. Rebolo, and S. Udry, "Statistical properties of 
exoplanets. II. Metalicity, orbital parameters, and space velocities," Astron. & Astrophys. 
398 363 (2003); 

C. Laws, G. Gonzalez, K. M. Walker, S. Tyagi, J. Dodsworth, K. Snider, and N. B. Suntzeff, 
"Parent Stars of Extrasolar Planets. VII. New Abundance Analyses of 30 Systems," As- 
trophys. J. 125 2664 (2003); 

G. Gonzalez, D. Brownlee, and P. Ward, "The Galactic Habitable Zone: Galactic Chemical 
Evolution," Icarus, 152, 185 (2001). 

M. Tegmark and M. J. Rees, "Why is the CMB fluctuation level 10" 5 ?," Astrophys. J. 



499, 526 (1998) |arXiv:astro-p h/9709058 1 . 

[29] M. Rees, J. Ostriker, "Cooling, dynamics and fragmentation of massive gas clouds: clues 
to the masses and radii of galaxies and clusters," Mon. Not. R. astr. Soc. 179 541 (1977). 

[30] K.Lang, "Astrophysical formulae, " Springer- Verlag (1974); 

M. Longair, "High Energy Astrophysics" Cambridge University Press (1981). 

[31] F. Hoyle, "On the fragmentation of gas clouds into galaxies and stars," Astrophys. J. 113 
513 (1953); 

M. Rees, "Opacity-limited hierarchical fragmentation and the masses of proto stars," Mon. 
Not. R. astr. Soc. 176 483 (1976). 

[32] M. Tegmark and M. Rees in [2E]; 

E. M. Leitch, and G. Vasisht, "Mass extinctions and the sun's encounter with spiral arms." 



30 



New Astronomy, 3 51 (1998) |astro-ph/9802l74| ; 

J. Garriga and A. Vilenkin, "Testable anthropic predictions for dark energy," Phys. Rev. 



D 67, 043503 (2003) |arXiv:astro-ph70 210358 ; 

M. Tegmark, A. Aguirre, M. Rees and F. Wilczek, "Dimensionless constants, cosmology 
and other dark matters," arXiv: astro-ph /05 1 1 774[ 

Ya. B. Zeldovich, "Gravitational Instability: An Approximate Theory for Large Density 
Perturbations," Astron. & Astrophys. 5 84 (1970); 

R. A. Sunyaev and Ya. B. Zeldovich, "Formation of Clusters of Galaxies; Protocluster 
Fragmentation and Intergalactic Gas Heating," Astron. & Astrophys. 20 189 (1972). 

T. Yanagida, in Proc. Workshop on the Unified Theory and the Baryon Number in the 
Universe, Tsukuba, 1979, eds. O. Sawada and S. Sugamoto, Report KEK-79-18 (1979); 
M. Gell-mann, P. Ramond and R. Slansky, in "Supergravity" North-Holland, Amsterdam, 
1979, eds. D.Z. Freedman and P. van Nieuwenhuizen. 

M. Fujii, K. Hamaguchi and T. Yanagida, "Predictions on the neutrinoless double beta 
decay from the leptogenesis via the LH(u) flat direction," Phys. Lett. B 538, 107 (2002) 
|arXiv:hep-ph/0203189| ; 

G. C. Branco, R. Gonzalez Felipe, F. R. Joaquim and T. Yanagida, "Removing ambiguities 
in the neutrino mass matrix," Phys. Lett. B 562, 265 (2003) |arXiv:hep-ph/0212 341 . 

M. Bolz, A. Brandenburg and W. Buchmuller, "Thermal production of gravitinos," Nucl. 
Phys. B 606, 518 (2001) |arXiv:hep-ph/0012052 . 



W. Buchmuller, K. Hamaguchi, M. Ratz and T. Yanagida, "Supergravity at colliders," 
Phys. Lett. B 588, 90 (2004) |arXiv:hep-ph/0402179|, "Gravitino and goldstino at collid- 



ers, 



arXiv:hep-p h/0403203 



M. Ibe, K. I. Izawa and T. Yanagida, "Realization of minimal supergravity," Phys. Rev. 
D 71, 035005 (2005) |arXiv:hep-ph/0409203| ; 

M. Ibe, T. Moroi and T. Yanagida, "Dark matter and baryon asymmetry of the universe 
in large-cutoff supergravity," Phys. Lett. B 620, 9 (2005) |arXiv:hep-ph/0502074| . 

M. Zaldarriaga, "Non-Gaussianities in models with a varying inflaton decay rate," Phys. 



Rev. D 69, 043508 (2004) |arXiv: astro-ph/ 0306006| ; 

D. Babich, P. Creminelli and M. Zaldarriaga, "The shape of non-Gaussianities," JCAP 
0408, 009 (2004) |arXiv:astro-ph/0405356 . 



M. Tegmark, A. Aguirre, M. Rees and F. Wilczek in 



31 



[41] V. A. Ambarzumian, Ann. Leningrad State U., No. 22 (1938); 

L. Spitzer Jr., "The stability of isolated clusters," Mon. Not. R. astron. Soc. 100 396 
(1940); 

S. Chandrasekhar, "Principles of stellar dynamics" University of Chicago Press (1942); 

L. Spitzer, Jr., "Distribution of galactic clusters," Ap. J. 127 17 (1958); 

L. Spitzer, Jr. and R. Harm, "Evaporation of stars from isolated clusters," Ap. J. 127 544 

(1958); 

I. King, "The escape of stars from clusters III," Astron. J. 63 465 (1958). 



32 



